arXiv: 1509.04798v2 [hep-ph] 30 Jan 2016 


Re+e- and an effective QCD charge 


J. D. Gomez^ and A. A. Natale^’^ 

^Centro de Ciencias Naturais e Humanas, Universidade Federal do ABC, 09210-170, Santo Andre - SP, Brazil 
^Instituto de Fisica Teorica, UNESP, Rua Dr. Bento T. Ferraz, 

271, Bloco II, OII 4 O-O 7 O, Sao Paulo - SP, Brazil 
(Dated: February 2, 2016) 

We consider the electron-positron annihilation process into hadrons R^+e- up to C1 (q?s) and we 
adopt the smearing method suggest by Poggio, Quinn and Weinberg to confront the experimental 
data with theory. As a theoretical model we use a QCD coupling constant frozen in the low energy 
regime, where this coupling can be parameterised in terms of an effective dynamical gluon mass 
(rrig) which is determined through Schwinger-Dyson equations. In order to find the best fit between 
experimental data and theory we perform a study, that, within the uncertainties of the approach, 
has a minimum value when nig / Aqcd is in the range 1.2 — 1.4. These values are in agreement 
with other phenomenological determinations of this ratio and lead to an infrared effective charge 
as(0) « 0.7. We comment how this effective charge may affect the global duality mass scale that 
indicates the frontier between perturbative and non-perturbative physics. 

PACS numbers: 12.38.Bx, 12.38.Aw, 12.38.Lg 


I. INTRODUCTION 

The standard model is one of the most successful the¬ 
ories of the last century, whose physical quantities are 
computed at the loop level with high precision showing 
remarkable agreement with the high energy experimental 
data. On the other hand, the theory that explains the 
strong interaction, known as Quantum Chromodynam¬ 
ics (QCD), is governed by an asymptotically free gauge 
field theory involving elementary quark and gluon helds 
at high energy or short distances. This behaviour at high 
energies is amenable to perturbative theory calculations 
using Feynman diagrams. This can be used to calcu¬ 
late physical mass-shell process. At low energy there is 
no justihcation for a perturbative treatment of QCD. In 
this regime non-perturbative approaches come into play. 

A non-perturbative approach for QCD is provided by 
the Schwinger-Dyson equations (SDE), whose study has 
revealed relevant progress in the recent years, most of 
them related to the disappearance of infrared (IR) di¬ 
vergences when the theory is resumed in a particular 
gauge invariant scheme, known as pinch technique [T]. 
The infrared QCD coupling turns out to be IR finite 
when gluons develop a dynamically generated mass (rUg) 
in this non-perturbative approach. This point was first 
demonstrated in Ref. [5] and was also discussed at length 
in Refs. [mni among many other references of this group. 
These results were checked independently by different lat¬ 
tice simulations [TTHI^ , has been studied in different ap¬ 
proaches [TM^ . and is, step by step, being accepted 
as a cure of the infrared QCD divergences. The hni- 
tude of the QCD coupling constant can be related to 
a non-perturbative IR fixed point, which is a property 
of dynamical mass generation in non-Abelian theories 
|23j . The phenomenological consequences of such IR fi¬ 
nite coupling, or non-perturbative fixed point, have been 
discussed in Ref. [2l] , and recently we have discussed how 
this non-perturbative fixed point can change the local 


minimum of a renormalisation group improved effective 
potential |25| . This change of minimum state may pro¬ 
duce noticeable modifications in the physical properties 
of the model studied in Ref. |26j . 

From the phenomenological point of view the theoreti¬ 
cal results leading to a finite QCD coupling is most than 
welcome, since many models describing strong interac¬ 
tion physics at low energy or small transferred momenta 
make use of an IR finite moderately small coupling con¬ 
stant. We present in the sequence a partial list of model 
calculations using an IR finite coupling constant: 1) The 
description of jet shapes observables requires an IR cou¬ 
pling equal to 0.63 [17], 2) Quarkonium potential models 
use an IR coupling of order 0.6 [H], 3) The theoretical 
ratio Re+e- can fit the experimental data with an IR cou¬ 
pling approximately 0.8 [H], 4) Calculation of quarko¬ 
nium fine structure in the framework of the background 
perturbation theory require a coupling as low as 0.4 m, 

5) QCD-inspired models describing total hadronic cross 
sections make use of an IR coupling of the order 0.5 EH, 

6) The experimental data on the unpolarised structure 
function of the proton is fitted with a coupling constant 
in the range 0.4 — 0.56 [32|. Besides the phenomenolog¬ 
ical indications, there are also several theoretical hints 
that the QCD coupling does not increase abruptly at low 
energies, and this fact may explain results claiming that 
the frontier between perturbative and non-perturbative 
physics may occur at relatively small momenta |33j . Of 
course, many of these calculations are based on models 
that cannot be fully connected to QCD, and even when 
related to QCD they may be, for example, renormalisa¬ 
tion scheme dependent, but they provide one hint about 
what we can expect for the IR value of the strong cou¬ 
pling constant. 

It should be remembered that a frozen coupling con¬ 
stant of the order that we described in the previous para¬ 
graph may contradict what is known about chiral sym¬ 
metry breaking in QCD 135] . It is indeed a problem 
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in the realm of Schwinger-Dyson equations to obtain the 
right chiral parameters for the value of the quark con¬ 
densate, /tt and other quantities when the gluon acquires 
a dynamically generated mass and the coupling freezes 
in the infrared. Possible solutions in this approach, in¬ 
voking non-perturbative QCD aspects, were pointed out 
only in the last years [3111371. The calculation of Ref. [51] 
is consistent with lattice data showing that the chiral 
and confinement transitions happens at the same tem¬ 
perature [38], and within this model we can obtain rea¬ 
sonable values for the chiral parameters [55] . Fortunately 
the problems raised by Peris and de Rafael [53] do not 
lead to a “no-go” theorem, and the solution of Ref. [51] is 
one possible way to evade these difficulties. 

In this work we will study the ratio Rg+e- '■ 


Re+e~ (®) 


a{e~^e —>■ hadrons) 
CT(e+e“ —>■ ’ 


( 1 ) 


as a tool to determine the infrared value of an effective 
QCD charge, which is related to the ratio nig/AqcD, 
where Aqcd is the QCD characteristic scale. A similar 
kind of analysis has already been performed in the known 
work of Mattingly and Stevenson [29] in the context of 
Optmised Perturbation Theory (OPT) and in Ref. [50] 
in the context of Analytic Perturbation Theory (APT). 
Here we will differ from them using an effective QCD 
charge obtained from Schwinger-Dyson equations, verify¬ 
ing that the experimental data is fitted only by a narrow 
range of mg/KqcD values. We recall that Rg+g-{s) can 
be calculated perturbatively, and can be developed as a 
power series in the QCD coupling, and in terms of the 
parameter 


a{s) = as(s)/7r. 


( 2 ) 


which, in our case, will be improved using the effective 
charge obtained through the SDE in an gauge invariant 
way using the Pinch technique [T]. 

The organisation of this work is the following: In Sec¬ 
tion II we discuss the hadronic cross section Rg+g-{s) 
in both partonic level and with massive quarks, and 
after that we use the non-perturbative effects in the 
QCD coupling constant to calculate the behaviour of ra¬ 
tio Rg+g-{s). In Section III the theoretical calculation 
is compared with the experimental data using smeared 
functions as proposed by Poggio, Quinn and Weinberg 
m, obtaining results that confirm the IR freezing of the 
QCD coupling constant. In Section IV we use the concept 
of global duality and discuss how this effective charge af¬ 
fects the scale that indicates the matching between the 
perturbative and non-perturbative physics. In Section V 
we draw our conclusions. 


II. HADRONIC CROSS SECTION A^+^-(s) 

The e“''e“ annihilation into hadrons is one of the most 
important processes for testing the theory of strong in¬ 
teraction. This process provides a fundamental QCD 


test, supplying evidence for the existence of colour [55] . 
At lowest order the total hadronic cross section is ob¬ 
tained by simply summing over all kinematically accessi¬ 
ble flavours and colours of quark-antiquark pairs. 


Rg+g- (s) — 


u{e^e —>■ hadrons) 
cr(e+e“ —^+/r“) 



(3) 


Here the qi denote the electric charges of the different 
flavours of quarks, and this is recognised as the parton 
model result. 

Real and virtual gluon corrections to this basic pro¬ 
cess ^ will generate higher-order contributions to the 
perturbative series. The second and higher-order correc¬ 
tions in perturbation theory were computed a long time 
ago in the zero quark mass limit in Ref. [53], and can be 
expressed order by order in a perturbative series 


"/ 

Rg+g-is)=3Y,<lfi^ + ns)), ( 4 ) 

i 

where R. has the form 

7?.(s) = a(l-I-rio-I-r 2 a^-I-• • •), (5) 

and depends upon a single kinematic variable s = , the 

c.m. energy. As shown in Eq. a is defined by the QCD 
coupling constant over tt. The coefficients of the Eq. ([^ 
depend of the renormalisation scheme. In the modified 
Minimal Subtraction scheme {MS) the coefficients were 
computed in Ref. [55] . 

ri = 1.986 - 0.1153n^ (6) 

r2 = -6.637 - 1.200n/ - 0.00518n^, (7) 

where n/ is the number of flavours. 

We will follow closely the pioneering work about OPT 
by Mattingly and Stevenson [55], but instead, we will 
consider the non-perturbative approach to the QCD cou¬ 
pling constant obtained through SDE. The approach de¬ 
scribed above is valid for massless quarks, nig = 0. To 
include quark mass effects we use the approximate result 

[551 El [55] 

Rg^g-{Q^) = iY.^f-^^^-vl)[l + g{vi)Tll (8) 

i 

with the sum over the active quark flavours that are 
above threshold, i. e. those with masses lower that Q/2 
{rrii < QI2), and 



= ( 
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Here Vi represents the quark velocity, so that Vi = 0 
corresponds to the heavy quark threshold, and the results 
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for massless quarks are recovered in the relativistic limit 

Vi —)■ 1 . 

For the theoretical calculation of i?e+e-(s) we shall 
take 7?.(s) giving by Eq. and ® with the coefficients 
ri and r 2 given by Eqs.([M and Q respectively, but as¬ 
suming an expansion in terms of the non-perturbative 
QCD coupling constant obtained from QCD Schwinger- 
Dyson equations in the pinch technique approach. A 
quite general expression for this coupling is given by [1] 


,{k^) = 


47r/3o In 




A2 

^^QCD 


-1 


( 10 ) 


where the function f {k"^^ m?{k'^)'j is determined as a fit 
to each specific value of the dynamically generated effec¬ 
tive gluon mass mg. Eor simplicity we have adopted the 
following expression for this coupling 




47 r/3o In 




-1 


( 11 ) 


where /3o = (lliV — 2n/)/487r^, Aqcd = A is the charac¬ 
teristic QCD scale, and the function mg{Q^) represents 
the dynamical gluon mass given by [45] 


mW) 


m„ 




( 12 ) 


The advantage of Eq. ( [11] ) is that we can vary the IR mg 
value without the neecim solving the coupled SDE for the 
gluon propagator, and as long as we consider the 3 < p < 
4 interval we accommodate the early Cornwall’s result 
[2|, the ones of Ref.[4| without noticeable differences, and 
this range is also consistent with the phenomenological 
values obtained in Ref.|31|. It is important to note that 
at high energies we recover the usual perturbative strong 
coupling. 

The effective QCD charge discussed in the previous 
paragraph has been obtained in one specific scheme 
(SDE and Pinch Technique) leading to a particular 
derivation of the non-perturbative QCD effective cou¬ 
pling. Although this scheme is gauge invariant [T] it may 
be claimed that we could have different definitions for 
the non-perturbative QCD coupling, and argue why we 
should consider Eq.(ll) as representative of the actual 


behaviour of the QCD charge. About this we can first say 
that Eq. 0 matches with the perturbative QCD cou¬ 
pling at high momenta. Secondly, lattice data and SDE 
solutions are clearly pointing to the existence of a dynam¬ 
ical mass scale for the gluon propagator, and if this is true 
we can prove the existence of a non-perturbative IR fixed 
point, i.e. an IR frozen conpling constant |23j . There¬ 
fore, Eq. 0 certainly matches the expected ultra-violet 
(UV) and IR behaviours of the QCD coupling constant 
no matter the scheme used to determine this effective 
charge, and we shall assume that in our calculation we 
can replace the perturbative coupling by this effective 


one. 


Going back to the R^+e- (s) ratio, one interesting issue 
is to compare the effects of quark masses at low energy on 
this ratio when we consider the effective QCD coupling 
constant as- In other words, we take both equations Eq. 
(j^ and Eq. (j^ at order 0{as), and compare them. So 
these equations becomes 




i?e.e-(Q^)=3 5:g?(l+“^) 


(13) 

(14) 


where as(fc^) is given by Eq. 0. We take standard 
values for the current-quark masses [IS]: = 2.4 MeV, 

md = 4.9 MeV, ms = 100 MeV, me = 1.3 GeV and 
A = 300 MeV. 



Q (GeV) 


FIG. 1. The one-loop theoretical ratio Re+ e- behaviour with 
the non-perturbative coupling of Eq.llll with p = 4 and 
mg IA = 1.2. The mass effects on the R^+e- is given by 


Eq.(|13[) (solid line), while the Eq. ( |14[ ) is represented by the 
dashed line. The calculation with the perturbative coupling 
constant is also shown by the dot-dashed line. 


From Fig.([^ we can see that the mass effect on the 
ratio Re+e- is negligible for light quarks, and become 
important with the opening of the charm threshold. Note 
that not only the freezing in the QGD conpling constant 
(Eq.(|Tl])) was considered into Eqs.(13) and (141, but the 
perturbative coupling constant has been considered into 
Eq. ( |I4| ) as well, and this one is represented by the dotted- 
dashed line in Eig.(|^. It seems that the mass effect, 
when we use the IR finite charge, is not so strong in the 
low energy regime (below 1 GeV) in comparison with the 
pure perturbative calculation. 

We have considered the expression for R{s) in the MS 
and the same happens for the quark masses. The mass 
effect is only relevant for heavy quarks, does not affect 
the main region of R(s) that we study when the coupling 
is infrared finite, and the mass effect was not included in 
the final result. However we do use a different scheme to 
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describe the infrared finite coupling! It has been argued 
that the Green’s functions obtained through the combi¬ 
nation of the pinch technique with the background field 
method (PT-BFM) are gauge invariant and renormaliza¬ 
tion group independent [37] (and so the coupling constant 
[48jb i.e. they are independent of any renormalization 
mass /i. This means that in principle we could obtain 
a coupling that would be independent of ambiguities in 
its determination. However this is not the case. The 
Schwinger-Dyson equation (SDE) for the gluon propaga¬ 
tor, from where it is obtained part of the information 
leading to the infrared coupling, has to be solved impos¬ 
ing that the non-perturbative propagator is equal to the 
perturbative one at some high-energy scale (/i), or com¬ 
paring the SDE propagator to the lattice data. After 
obtaining the QCD propagators we can determine the /i 
independent coupling through one specific relation of two 
point correlators. This procedure is not unambiguous 
and contains all the numerical uncertainty related to this 
specific calculation. We may argue that the mass genera¬ 
tion mechanism in the PT-BFM approach minimizes the 
vacuum energy P| and should reflect a scheme indepen¬ 
dent quantity, however the full procedure to obtain the 
infrared finite coupling contains truncation and numeri¬ 
cal approximations that are hard to be estimated at the 
present level of the SDE solutions. As a matter of com¬ 
pleteness we recall that another possible determination of 
the IR-behavior of as is based on Light-Eront holographic 
QCD [49|. The only assumption that underlies Light- 
Front holographic QCD is that QCD is conformal in the 
IR, in a procedure claimed to be renormalisation scheme 
independent m- Similarly to the above discussion, the 
determination of the coupling uses a particular definition, 
which is called the gi scheme, that may introduce uncer¬ 
tainties in the same way as the PT-BFM approach. The 
gi scheme can be related to other couplings, as the one 
assumed in our work, using fundamental QCD relations, 
and these relations suggest that, once transformed to the 
one discussed here, the coupling freezes near 0.7 m, i-e. 
the same value found by us. 


III. SMEARING Ae+^-(s) 

The comparison of the theoretical prediction with 
the experimental data for Re+e- is not possible rig¬ 
orously, because there is no direct correspondence be¬ 
tween the perturbative quark-antiquark thresholds and 
the hadronic thresholds and resonances of the data. The 
only way to do this is using the smearing method pro¬ 
posed by Poggio, Quinn and Weinberg mi- Although 
this method was thought to be used at high energies, in 
the Refs. jUj and HD] it was used at low energies in 
different contexts and with interesting results. 

The experimental data that we are going to use has 
been taken from Particle Data Group R ef. [46j. Of par¬ 
ticular interest to us is the region from Q = 0 up to 6 
GeV. The data and errors are shown in Fig.([^. In that 


figure we have done a zoom to better show the resonances 
in the region of interest to us. 
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FIG. 2. World data to ratio R^+^-{s) = a(e~^e~ —>■ 
hadrons, s)/(T(e"'"e“ —>■ a{e^e~ —>■ hadrons, s) 

is the experimental cross section corrected for initial state 
radiation and electron-positron vertex loops, a{e^e~ —>■ 
,s) = 47 ra^(s)/ 3 s. Data errors are total below 2 GeV 
and statistical above 2 GeV. 



Figure ^ shows our data compilation, up to 6 GeV. 
The red solid line represents our fit, where the quite nar¬ 
row resonances p, </>, J/i^, V'(3686), 1/1(3770), as assumed 
in Ref. |29] . were not included into the data compilation. 
The data go well beyond the b quark threshold, but they 
have no real effect on the results that we shall present. 
In the fit of the Fig. (§ are included four resonances, 
the oj, p', i/i(4040), i/)(4415). The red curve was obtained 
with the following fit: 


E 


(Q2 _ Mfy + A, 


-l- C -l- DQ^, 


(15) 


where Mi = 0.781, M 2 = 1.65, M 3 = 4.04, M 4 = 4.42 
GeV^ and we summarise the values and errors of the 
fit parameters in table (|T]). To perform the fit we used 



Estimate Standard Error 

Ai 

3.61 X 10"^ 9.80 X 10“^ 

A 2 

4.50 X 10"^ 6.73 X 10“® 

^3 

3.42 X 10"® 5.69 x 10“'‘ 

A 4 

1 . 11 x 10 "^ 2.26x10“® 

Bi 

5.37 X 10“® 3.55 x lO”’’ 

B 2 

9.11 X 10“® 4.80 X 10“'^ 

B 3 

1.20x10“® 4.87x10“’’ 

Bi 

8.71 X 10“® 4.95 X 10“^ 

C 

-4.63 X 10“® 7.00 X 10“^* 

D 

6.66 X 10“® 2.26 x 10“® 


TABLE I. Fit parameters and errors respectively leading to 
the fitted curve in Fig. li¬ 
the NonlinearModelFit package of Mathematica software 
obtaining a value for the number of adjusted points 
equal to R? = 0.9921. A quite detailed discussion about 
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the errors appearing in the analysis of the experimental 
data is already presented in the second work of Ref. [29] . 



Q (GeV) 


FIG. 3. Compilation of experimental i?e+e- data (excluding 
narrow resonances). The red solid line represents the fit of 
Eq.pl. 


In order to compute the smeared experimental data 
we employ the same approach as the one of Ref. . We 
assume that the narrow resonances have a Breit-Wigner 
form 


RrcsiQ^) 


B 

y2 « '‘(Q2_^2)2+JVf2r2’ 


(16) 


where a is the QED coupling, and M, F, Ba, Bh are 
the mass, width, lepton branching ratio, and hadronic 
branching ratio respectively. The narrow resonances that 
were excluded from the fit can now be represented by a 
delta function through the follow transformation 


(Q2_M2)2 + M2r2 




(17) 


For these we can write the Eq.(16) as 


RresiQ^) « ^^^BuBhSiQ^ - M^). (18) 


In order, to compare both the theoretical and exper¬ 
imental data we used the “smearing” procedure applied 
by Poggio, Quinn, and Weinberg (PQW) in Ref.|4T|. 
There they defined the “smeared” ratio by 


Rpqw{Q^-,^) ^ _q2^^2 ^^2- ( 19 ) 

The smeared ratio could be written as m 

2iRpQw{Q^; A) = n(g 2 + ^a) - b[q'^ - zA). ( 20 ) 

where n(z) is the vacuum-polarization amplitude. The 
best choice of A is the smallest value that will smooth 
out any rapid variations in either the experimental or the 
theoretical i?e+e-- It turns out that this depends upon 
the energy region one is interested in. Around charm 


threshold a A = 3 GeV^ or more is necessary, while in 
the lowest-energy region a A as small as 1 GeV^ can be 
used mj. The idea now is to apply this smearing to 
both the theoretical and experimental i?g+e-’s and then 
compare the results. 

To determine the theoretical smeared ratio we 


can solve the Eq.(19l by numerical integration. In order 


to deal with the different thresholds we performed the 
numerical integration considering several intervals, from 
ruu to Qmax and over the range 0 to mu we considered 
Re+e- = 0- We assumed Qmax = 6 GeV. From g„iax 
to oo, Re+e- remained constant, and in this region we 
can make the integration of Eq.(19) analytically. The 


quite narrow resonances are given by Eq.(18), therefore 


the Eq.( 19) can be solved without trouble, and their con¬ 


tributions become 


Rr 


9BuBhAMT 

a2[(Q2_M2)2-pA2]' 


( 21 ) 


The other resonances (w, p', ■!/’(4040), '0(4415)) were con¬ 
sidered in the fit with a Breit-Wigner function. 

The fitted experimental data was integrated numeri¬ 
cally as shown in Eq.([T^. We computed the “smeared” 
quantity for four different values A =1 GeV^, 1.5 GeV^, 2 
GeV^, and 3 GeV^. For simplicity we present in Fig.Q 
the result only for A = 1.5 GeV^. The smeared theo¬ 
retical Re+e- was computed, threshold by threshold, for 
values of mg/A between 0.7 and 2.4 and p = 4, and the 
values are also shown in Fig.Q. Note that the QGD 
coupling constant Eq.(ll) depends on mg/A, and conse¬ 
quently the ratio Re+e- has the same dependence. The 
shaded area in the figure was determined assuming ±7 % 
uncertainty around the experimental central value, ob¬ 
tained when we considered ad hoc variations of the many 
parameters in our fitting procedure and the possible nor¬ 
malisation errors pointed out in B.ef. |29j . 

In the Fig. most of the theoretical smeared lines 
are within the shaded area except for the result obtained 
with mg/A = 0.7 (magenta line with square data points), 
which is out of the shaded region. This fact will hap¬ 
pens for any mg/A value smaller than 0.7, which lead to 
large values of the coupling constant where the expansion 
in Re+e- (s) does not make sense anymore. This mean 
that there is not agreement between data and theory for 
mg/A = 0.7 and lower values within the assumed un¬ 
certainty. The theoretically smeared Re+e- with mg/A 
between 0.8 and 2.4 are within the shaded region. How¬ 
ever, larger mg/A values lead to theoretical Re+e- curves 
increasingly away from the experimental data, but with 
lines that barely can be distinguished from each other 
and we did not include these results in the figure. 

We have applied the smearing method for the ratio 
Re+e- with p = 4 and different values of A in order to 
compare experiment and theory. In principle, following 
Ref. |29| . a smaller A value is better to describe the low 
energy limit, and we shall limit ourselves to Re+e- up to 
1.5 GeV^ as shown in Fig. In order to compare the 
data we performed a test for our physical quantities. 
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(GeV) 


FIG. 4. Smeared R^+e- for ^ — 1-5 GeV^. The experimen¬ 
tal result is shown in the middle of the shaded region. We 
have drawn this shaded area assuming an error of 7% above 
and below the smeared experimental data, what helped us 
to determine different values of the ratio mg/A that are near 
this region. For example, the theoretical smeared Rf,+e- for 
mg/A = 0.7 is shown by the magenta line, which is out of the 
shaded area. 


provides the best match with the experiment^ and points 
out for one specific infrared value of the QCD coupling 
constant. It is important to remember that in the SDE 
equations mg is one input parameter and lattice data has 
not enough precision to pinpoint this mass scale, reveal¬ 
ing the interest on phenomenological determinations of 
this quantity. 



We choose a set of of measured observables and 

compute {Ol^{9)}, where 9 are parameters. Then we 
minimise the function 


xHS) = E 

i 


(or-o!'^(9)f 

(A0-p)2 ’ 


where the observable of our interest is the ratio Rg+e-, 
and the parameter is the dynamic gluon mass mg/A. 
With 21 points (degrees of freedom) and a deviation of 
data about 10% we develop a test for comparison of 
the smeared Rg+e- at ©(of). In the Table (|n| we set 
p = 4 and write the x^ values for several mg/A. The val¬ 
ues of rm/A are in the interval [0.8, 2.4], and according 
to Fig.Q they fall into the shaded region. 



A = 1.5 A = 2 A = 3 

00 

o 

II 

< 

1.37763 1.85683 2.70305 

mg/A = 0.9 

0.837894 1.24134 2.10877 

mg/A = 1.0 

0.712088 1.06814 1.92454 

mg/A — 1.2 

0.657517 0.98275 1.82944 

mg/A = 1.4 

0.666606 0.99339 1.84067 

mg/A — 1.6 

0.698763 1.03745 1.88845 

mg/A = 1.8 

0.744845 1.09813 1.95303 

mg/A — 2.0 

0.801059 1.16891 2.02676 

mg/A = 2.2 

0.865123 1.24620 2.10588 

mg/A = 2.4 

0.935314 1.32809 2.18823 


FIG. 5. x^ for the smeared quantity Rg+g- {Q^, A) with p = 4 
and different A values. 


As we discussed after Eq.(10l for each value of m^/A 
we have a slightly different SDE solution, whose differ¬ 
ences can be parameterised in terms of different p val¬ 
ues. Therefore we have changed the A values, the p val¬ 
ues and studied the x^ distribution. With p = 4 and 
A = 1.5, 2, 3 GeV^ we have found that the minimum is 
located at mg/A = 1.2. However, with A = 1.5 GeV^, 
A = 2 GeV^ and taking p = 3,3.5,4, we can see in 
Fig.([^ and Fig.Q that the minimum points are now lo¬ 
cated between mg/A = 1.2 and mg/A = 1.4. Note that 
the minimum for A = (1.5, 2) GeV^ and p = 3, 3.5 are 
approximately at mg/A = 1.4. Unless some SDE solu¬ 
tion deviates grossly from the results that are presently 
found in the literature, mg/A values between 1.2 and 1.4 
and p values between 3 and 4 lead to an infrared value 
of the coupling constant, according to Eq.(|TI|) with two 
quark flavors, approximately of 0(0.7). The minimum 
X^ seems to be dependent on the product pm^/A. 


IV. GLOBAL DUALITY AND THE INFRARED 
QCD COUPLING 


TABLE II. x^-values for Rg+g- with p — 4 and different A 
values (in units of GeV^). 

One interesting point in the table 0 is that for dif¬ 
ferent mg/A values and p = 4 all minimum values of x^ 
occur for m^A — 1.2. This behaviour can be better 
seen in Fig.(™, where v is the number of degrees of free¬ 
dom. This means that there is one value of mg/A that 


It is usually assumed that there is a frontier between 
perturbative and non-perturbative QCD at one specific 
mass scale, and there are proposals about how this scale, 
that separates long and short distances, can be deter¬ 
mined. One of these proposals is established by the con¬ 
cept of global duality and is explained in Refs.[52l l53] . 
The determination of this scale (sq) consist in doing a 
matching between perturbative QCD (pQCD) and ex- 
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4.5 



°'B.5 1.0 1.5 2.0 2.5 3.0 

rtig! h. 


FIG. 6. for fhs smeared quantity R^+f,- {Q^, A) with A = 
1.5 GeV^ and different p values. 



FIG. 7. x^ for the smeared quantity R^+f,- (Q^, A) with A = 
2 GeV^ and different p values. 


periment by the expression 


rso I rso 

/ dt-ImIV{t)exp= / dt-ImIV{t)pQCD- (22) 

Jo tr Jo ’r 

The So scale was determined to be of sq « 1.5 GeV^ 


Eq. (22) also imply 


ds R^xp / dsR^fi^Qj-^ 


(23) 


and we can use our previous result to determine sq. It 
is clear that we shall have a different sq value for the 


matching prescribed by Eq. (22) if we use the effective 


charge that we discussed before. This happens because 
in the case of perturbative QCD the right-hand side of 
Eq.(22| starts growing fast as we approach the Landau 


pole of the perturbative coupling and the solution ap¬ 
pears at a large mass scale. On the other hand assuming 
global duality for the IR finite effective QCD charge we 
do not expect a fast rise of the coupling at low momenta 
and consequently a slow increase of the right-hand side 
of Eq. (231 as we go to small momenta. 


To determine the scale sq we can define the quantity: 


e(Q ; Sq) — / ds (^Rexp Rtheor)- (24) 

JO 

and look for zeros of this equation. In the above equation 
we have a Sq and to^/A dependence besides the one on 
A and p. Eor simplicity we can solve Eq. (24) in the limit 
—>■ 0, for specific values of the smearing factor A, p 


and rag/K and look for a zero of Eq.(24). 


In the case of A = 1.5 GeV^, p = 4 and m-g/A = 1.2 we 
show in Table (fill that e((5^ —)■ 0) « 0 when Sg = 0.87 
GeV^. This, as discussed previously, is the result that 
we were expecting when we change the perturbative cou¬ 
pling by the effective charge of Eq. 0 , showing that 
the so called frontier between the perturbative and non- 
perturbative physics occurs at one smaller mass scale 
than the one determined with perturbative QCD. We also 
show in Eig. § a plot of e{Q^ —>■ 0) as a function of sg 
and rrig/A in the case of A = 1.5 GeV^ and p = 4, where 
we can show that the smallest e(0) values are obtained 
at one sg smaller than 1.5 GeV^. Actually the effective 
coupling that we are using has non-perturbative informa¬ 
tion, therefore we should not be allowed to say that sg 
is the mass scale indicating the transition between short 
and large distances, but just say that if we improve the 
QCD calculations with such effective charge we proba¬ 
bly could perform “improved perturbative” calculations 
a little bit deeper into the IR region. 


So 

e(0) 

0.6 

-0.0987 

0.8 

-0.0276 

0.87 

-0.0060 

1.0 

-0.0291 

1.2 

-0.0812 

1.4 

-0.1264 

1.6 

-0.1545 


TABLE III. e(0) with rUg/A = 1.2, p = 4 and different so 
values. 



2.0 


^IG. 8. e{Q^ —> 0) surface (Eq.(24l) calculated up to 0{al 
vith A = 1.5 and p = 4. 
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V. CONCLUSIONS 

We have studied the electron-positron annihilation 
process into hadrons i?e+e- up to C’(af), where the QCD 
coupling constant was described by an effective charge 
obtained in solutions of the Schwinger-Dyson equations 
in the pinch technique scheme. This effective charge is 
frozen in the infrared region and the frozen value is re¬ 
lated to a dynamically generated gluon mass. The exis¬ 
tence of such dynamical gluon mass, or an effective mass 
scale for the gluon propagator, has not only been ob¬ 
served in SDE solutions but has been confirmed by lat¬ 
tice simulations, and in this case we can prove that the 
coupling constant freezes in the IR region |23j . The main 
purpose of the work was the determination of the in¬ 
frared value of the dynamical gluon mass, related to the 
IR value of the effective charge, since this value enters 
as an input into the numerical solutions of the SDE, and 
lattice data including dynamical quarks can only obtain 
rough approximations to this quantity. 

To compare the Re+e- experimental data to the theo¬ 
retical calculation we adopted the smearing method sug¬ 
gest by Poggio, Quinn and Weinberg. In order to find the 
best fit between experimental data and theory, we per¬ 
formed a study of the resulting curves, that, within 
the uncertainties of the approach, leads to a minimum 
value when rngjAgcD is in the range 1.2 — 1.4. These 
values are in agreement with other phenomenological de¬ 
terminations of this ratio and imply an infrared effec¬ 
tive charge 0^(0) « 0.7, what is also in agreement with 
the result obtained by Mattingly and Stevenson when 
analysing Rg+g- in a different scheme [^ . 

The study indicates that the result is stable and 
independent of the smearing parameter A. All uncer¬ 
tainties about the full procedure have been discussed in 
Ref. [in]: and the main novelty is the use of the effective 
charge related to the dynamical gluon mass scale. Al¬ 
though there are many evidences for an IR finite gluon 
propagator and coupling constant, there is not an unique 
definition of the non-perturbative QCD charge. However, 
the effective coupling that we have considered here can 
map any possible behaviour of the QCD charge, i.e. is 
compatible with the UV behaviour predicted by pertur¬ 
bative QCD and freezes in the IR region as should be 
expected when the theory develops a dynamically gener¬ 
ated mass. 

Taking advantage of our Re+e- calculation we have dis¬ 
cussed what happens with the scale determined by the 
concept of global duality jSH [53] when the perturbative 
coupling is exchanged by the effective charge discussed up 
to now. We verify that sq, the mass scale indicating the 
transition between short and large distances, is changed 


to a smaller value. However, as should be understood 
from our procedure, we are assuming an “improved” cou¬ 
pling, and this fact has clear phenomenological conse¬ 
quences. One is that we could perform “improved per¬ 
turbative” calculations possibly going deeper into the IR 
region, although the main point is that it is necessary to 
build a bridge between perturbative QCD and the SDE 
results, which not only show a cure for the QCD IR di¬ 
vergences, but seems to provide a soft transition between 
the perturbative and non-perturbative QCD regimes. 

Although the idea of dynamical gluon mass generation 
was put forward a long time ago, leading to infrared finite 
Green’s functions |2] , due to the involved field theoretical 
aspects of the combination of the pinch technique with 
the background field method necessary to obtain gauge 
invariant and renormalization group independent quan¬ 
tities, only after the appearance of many lattice QCD 
simulations this subject had a revival, in such a way 
that very recent theoretical papers are revisiting the de¬ 
termination of a finite gluon propagator and coupling 
constant in three |54| and four dimensional QCD [47], 
trying to spread the ideas of this mechanism. On the 
other hand this subject also raised attempts to describe 
these infrared finite quantities in terms of an effective 
Lagrangian valid for infrared QCD (see Refs. [551157] and 
references therein). These attempts, although introduc¬ 
ing a hard gluon mass and leading to a “soft” BRST sym¬ 
metry breaking, seem to be compatible (at the “pertur¬ 
bative” level) with many lattice data about the freezing 
of the infrared QCD quantities. Therefore, including the 
discussion presented in this work, there are several new 
[3HHHT] and old [32] tests for the phenomenological con¬ 
sequences of an infrared finite QCD coupling and gluon 
propagator, and, as the time goes on, all these attempts 
may provide a new consistent view of infrared QCD. It 
is also interesting to note that in a recent work [63] it 
was determined the transition scale between perturbative 
and non-perturbative QCD using Hadron-Parton duality, 
which, when translated to the MS scheme, would give a 
value close to the coupling determined by us. This can 
only indicate that the different approaches lead to a re¬ 
sult that is quite stable. 
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